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Abstract 

Dynamical fermions induce via the fermion determinant a gauge-invariant 
effective action. In principle, this effective action can be added to the usual 



gauge action in simulations, reproducing the effects of closed fermion loops. 
^ ■ Using lattice perturbation theory at finite temperature, we compute for stag- 

o ' 

■ gered fermions the one-loop fermionic corrections to the spatial and temporal 

On _ 

plaquette couplings as well as the leading Zfq symmetry breaking coupling. 
A. Hasenfratz and T. DeGrand have shown that j3 c for dynamical staggered 

■ fermions can be accurately estimated by the formula j3 c = /?P ure — Af3p where 

■ A/3j? is the shift induced by the fermions at zero temperature. Numerical and 
analytical results indicate that the finite temperature corrections to the zero- 
temperature calculation of A. Hasenfratz and T. DeGrand are small for small 
values of k = t^Ji hut become significant for intermediate values of k. The 
effect of these finite temperature corrections is to ruin the agreement of the 
Hasenfratz-DeGrand calculation with Monte Carlo data. We argue, however, 
that the finite temperature corrections are suppressed nonperturbatively at 
low temperatures, resolving this apparent disagreement. The Zjy symmetry 
breaking coupling is small; we argue that it changes the order of the transition 
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while having little effect on the critical value of (3. 
PACS number(s): 12.38.G, 12.38.Mh, ll.10.Wx 



I. INTRODUCTION 



It has been known for some time that the effects of heavy dynamical fermions can be 
included in Monte Carlo simulations by a hopping parameter expansion of the fermion 
determinant. This is reminiscent of the Euler-Heisenberg Lagrangian of perturbative QED, 
in which the effects of electron loops are included in a gauge-invariant effective Lagrangian. 
Recently Hasenfratz and DeGrand have performed a zero-temperature calculation of the 
shift in the lattice gauge coupling constant /3, defined as (3 = 2N c /g 2 , induced by staggered 
dynamical fermions and applied the result to the finite temperature phase transition in 
QCD. Their result for the shift in the critical coupling, in the form f} c = /3,F ure — A/3p, was 
found to hold rather well down to small values of the fermion mass. It is convenient to work 
with a hopping parameter k defined by k = 2^;, where mj? is the fermion mass. As shown 
in Fig. [j], the Hasenfratz-DeGrand results for N t = 4 are in excellent agreement out to at 
least k = 2, and in reasonable agreement at k = 5; this is particularly surprising since Aflp 
is calculated using lattice perturbation theory at zero temperature. In order to understand 
the effects of finite temperature, we have calculated the one-loop fermionic corrections to 
the spatial and temporal plaquette couplings, as well as the leading Zn symmetry breaking 
coupling. 



II. RENORMALIZATION OF AT FINITE TEMPERATURE 

A. Perturbation Theory for A/3 

The 0(A 2 ) term in the gauge field lattice action including the one- loop finite temperature 
fermionic correction is given by 



S, 



off 



-9 2 E ^ j_] ^^Mp)Mp)4r[ D{ >) + D ^p) - D $] } ( 21 ) 



where 



D$(P) = 4[^E*in 2 (f ) - sin(f ) sin(^ 



(2.2) 
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Di » = \^W t LI f^W^ + f )^K^ + P f)S-\k + p)) , (2.3) 



and 



^-^E^/J^Iow^*)] (2.4) 

with the vertex functions given by 

= ijf, cos(A; M ) (2.5) 

and 

Q(*7«) = sin(fc M ) (2.6) 

with no sum over //. The inverse fermion propagator is 

S{k) = ^- + ij fM sin{k tl ). (2.7) 

This formula is a straightforward consequence of the lattice Feynman rules, which are given 
in Fig. |2|. The diagrams contributing to the fermionic renormalization of A/3 are shown in 
Fig. |3|. The first diagram, corresponding to D^°\ is the free lattice gluon propagator. The 
second diagram, corresponding to involves R and S only and survives in the continuum 
limit; note that R is the lattice analog of the continuum gluon vertex. The third diagram, 
corresponding to D^ 2 \ is a lattice tadpole diagram, and involves the vertex function Q, which 
is a feature only of the lattice theory. At zero temperature, this tadpole contribution vanishes 
after integration by parts [J3J. Finite temperature enters into the calculation only through 
the replacement of the integration over the ko variable appropriate for zero temperature by 
the sum over Matsubara frequencies 

27m 

ko = ^r (2.8) 
where n is integer-valued and T is the temperature. 
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B. Ward Identity at Finite Temperature 



At finite temperature, the term in Eq. Q2.1J) is necessary in order to show that the 
lattice form of the Ward identity 



'2V 







(2.9) 



still holds at finite temperature, even though the four dimensional hypercubic symmetry is 
broken. To show this, we first note the two identities 



S(k +p)- S(k) = 2 sin(^) R(k, + ^ 



and 



R(k, + m- R(k, - f ) = 2 sin(^)Q(U 



Use of the first identity gives 



1 



1 



<P^ 



[r{K + ^)S-\K+p,) - R(k, + ^S-'ik,)} }, (2.12) 



(2.10) 



(2.11) 



and after a simple shift of variables, use of the second identity yields the desired cancellation. 



C. Finite Temperature Decomposition of the Propagators 

At finite temperature there are two independent symmetric tensors of order p 2 which are 
four-dimensionally transverse In what follows, all expressions will be in the thermal rest 
frame. The first of the corresponding lattice tensors is specified by 

P!$(p) = P&\p) = (2.13) 

P s 

and the second is 
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P!$(P)-P$(P) (2-15) 

where 

P^( P ) = ^- P -f i (2-16) 
The lattice quantities p are defined by 

^ = 2sin(|) (2.17) 

P 2 = T,PI (2-18) 
/* 

p1 = £p1 (2-19) 

i 

The existence of these two independent tensors leads to separate renormalizations of 

the spatial and temporal gauge couplings at finite temperature. The first tensor, Pj£)(p), 

is associated with the magnetic, or spatial, part of the action, while the second tensor, 
Puv(p) ~ Pju)(p)i i s associated with the electric, or temporal, part of the action. 

We find that 

d 3 k 



and 



where 



1 r 77 n h 

Aft ^^-/^ft 1,2) P.20) 
AA-^E^/JII^M.I) (2.21) 



$(Jfe; [I-, v) = 32B- 2 (k) cos 2 (£; M ) cos 2 ^) 

-4096£ _4 (fc) sin 2 (^) cos 2 (A; M ) sin 2 (^) cos 2 (k u ) (2.22) 



and 



5(A;) = ^ + 4^sin 2 (A; a ). (2.23) 

K a 

As the temperature is taken to zero, the two expressions smoothly approach each other to 
give the zero-temperature result. 
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D. Numerical Results for A/3 



The integrals ( 2.20 ) and ( |2.21| ) were evaluated numerically by calculating mode sums 



for large values of N s and various values of N t . Figure |] compares A(3 per fermion (spatial 
and temporal) vs. k for N t = 4 and N t = 6 with the zero temperature result of Hasenfratz 
and DeGrand JTJ . The finite temperature values approach the zero temperature result from 
below as a consequence of the antiperiodic boundary conditions; as N t increases, the sum 
includes more terms in the region near p = 0, which dominates for small mp. As might be 
expected, the temporal shift in j3 is more sensitive to the effect of finite temperature than 
the spatial shift. For small values of k, corresponding to large values of the fermion mass, 
the effects of finite temperature are small. This is easily understood, since k is much smaller 
than N t . However, for intermediate values of k, finite temperature corrections ruin the 
excellent agreement between the zero-temperature calculation and the Monte Carlo results 
discussed in Sec. Q. 

E. Image Expansion 

The connection between the the zero and finite temperature result can be understood 
more physically by transforming the sum over Matsubara frequencies into a sum over images 
using the Poisson summation formula for antiperiodic boundary conditions 

E w F ^ = B- 1 )" / * ^r F (poy nNtpo ( 2 - 24 ) 

so that, for example, the shift in the spatial coupling is given by 

r * d 4 k 00 r * d^h 

Ap s = -N c ^ k ;l,2)-2N c J2(-l) n Z *(*; 1, 2) cos(nN t k ) (2.25) 

J-n (2tt) 4 n=1 J-n (2tt) 4 

with a similar result for Aj3 t . This form has a simple physical interpretation: the first 
integral is the zero-temperature shift, and the integer n in the second term labels the net 
number of times the fermion wraps around the lattice in the temporal direction. The finite 
temperature corrections result from the 0(A 2 ) expansion of image diagrams such as those 



depicted in Fig. |5|. Numerically, the dominant corrections to the zero-temperature result 
come from the first few values of n, with the n — 1 and n = 2 terms accounting for more 
than 90% of the finite temperature correction for k < 2.0. 

Although not apparent in our perturbative calculations, in order to maintain gauge 
invariance, the vertical segments of the image diagrams must be accompanied by powers of 
Polyakov loops. It is an observed feature of simulations with dynamical fermions that the 
symmetry is approximately maintained in the low-temperature phase. This suggests that 
the image contibutions may be negligible below f3 c . Thus, the zero temperature corrections 
to (3 are suppressed nonperturbatively in the confined regime. In particular, just below f3 c 
the zero-temperature result will hold for A/3. Figure ^| illustrates an idealized behaviour for 

^/^fermion ^ 

function of /3 pure . 



III. Z N SYMMETRY BREAKING IN THE EFFECTIVE ACTION 

There is another set of terms induced by the fermion determinant only at finite temper- 
ature. As is well known, the symmetry of the pure gauge theory is explicitly broken 
by dynamical fermions. To lowest order in the hopping parameter expansion, a path of 
N t hops around the lattice in the temporal direction produces a effective coupling to the 
Polyakov loop, explicitly breaking the Zn symmetry. Evaluating the fermion determinant 
in a constant A$ background field and applying the Poisson summation formula [Eq. ( [2.241) 1 
once again, we find an additional contribution to the effective action: 



S, 



cff 



EEK 

x n=l 



n+l 



h n (n)Re 



TrP n (x) 



(3.1) 



where P(x) denotes a Polyakov loop and the couplings h n are given by 



-4N t 



■ A 

-7T <xF 



In 



Ak 2 



+ I>in 2 (g At ) 



cos(nN t q ) 



(3.2) 



This result can also be obtained by using a contour integral to evaluate the sum over Mat- 
subara frequencies. The leading term in this effective action has been discussed for the case 
N t = 2 §. 
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A. Numerical Results for h 



The maximum values of the h n are obtained when mp = 0. For N t = 4 we find /i™ ax = 
0.107, /i™ ax = 0.00445, and increasingly smaller magnitudes for higher order couplings. 
Because h± favors Zjq breaking, it acts to lower the critical value of (3. Unlike the D^ l \p) 
term considered in Sec. [TT], this effect cannot be directly included as a finite shift in (3. The 
most direct way to determine the shift in (3 due to hi term is to perform a Monte Carlo 
simulation of the pure gauge theory with an hi term added to the action. The numerical 
simulation results presented in this section were obtained from runs of 40,000 sweeps (after 
thermalization) on a 10 3 x 4 lattice using a variety of work stations. 

We have observed that this additional source of A/3 is small compared to the renormal- 
ization of the plaquette couplings discussed in the preceeding section, regardless of m F . For 
example, even hi — 0.1 at N t — 4 leads to a shift in (3 per fermion of 0.00325. This value 
of hi corresponds to mp = 0.17 which yields a zero-temperature predicted shift in (3 per 
fermion of 0.104. Although the effect of the hi term in the effective action on the critical 
value of (3 is quite small, hi has a profound effect on the character of the transition. Figure [7] 
shows the frequency distribution for the spatial expectation value of the Polyakov loop as 
a function of hi at the appropriate (3 c {hi) on a 10 3 x 4 lattice. As hi increases, the peaks 
associated with the two phases move closer together until they merge at what is presumably 
a second-order critical end-point. A first-order phase transition exists for values of hi smaller 
than approximately 0.08. 

The Monte Carlo results for the phase diagram of a pure gauge theory with an hi 
coupling can be mapped back to the phase diagram of QCD with dynamical quarks by 
including both the hi term and the shift A/3. Including zero-temperature plaquette coupling 
renormalization, the endpoint of this first-order phase transition line maps to the point 
(0.394, 4.68) in the (mp, /3)-plane for the case of sixteen degenerate staggered fermion species. 
This is roughly consistent with the endpoint of the first-order phase transition observed in 
simulations with dynamical staggered fermions 0. However, a definitive comparison will 
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likely require significant computational resources. 



IV. CONCLUSIONS 

Figure |8|, a graph of j3 c versus k, summarizes the results we have obtained. As shown 
by Hasenfratz and DeGrand, the zero-temperature shift in the coupling constant due to 
dynamical fermions nicely accounts for the shift in j3 c . We have found that finite temperature 
corrections to the gauge coupling renormalization lift the degeneracy of the spatial and 
temporal couplings, and the results are significantly different from the zero-temperature 
results. As can be seen from the figure, they are in conflict with the Monte Carlo data. 
At finite temperature, dynamical fermions couple to Polyakov loops via loops circling the 
lattice in the timelike direction. This also acts to shift the value of j3 c . This shift is small, 
however, and does not restore the success of the zero-temperature calculation. 

As we have shown, the image expansion makes it plausible that the success of the zero- 
temperature calculation in determining the critical value of (3 is due to a nonperturbative 
suppression of finite temperature effects in the low temperature regime. Specifically, the 
small expectation value of the Polyakov loops at low temperatures indicates a suppression 
of those quark paths which account for finite temperature corrections. This leads to the 
behavior of the effective gauge coupling constants shown in Fig. |5] above. If this picture 
is correct, a gauge theory with dynamical fermions of intermediate mass is well modeled 
by a pure gauge theory in which the spatial and temporal couplings are slightly different. 
Along the critical line, both couplings jump discontinuously. This could be studied in more 
detail by adjusting spatial and temporal couplings in a pure gauge theory simulation so that 
plaquette expectation values matched those observed in dynamical simulations. This might 
also afford an opportunity to consider the effects of finite spatial sizes. Such effects are easily 
included in Eq. p.25| for the coupling shift by choosing the spatial mode sums appropriately. 



While the coupling to Polyakov loops induced by dynamic fermion loops seems to play 
little role in determining (3 C , this coupling does influence the order of the transition; in our 
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N t — 4 simulations of a pure gauge theory with an additional coupling hi to Polyakov loops, 
a sufficiently large value of hi causes the first order line to terminate, while shifting f3 c very 
little. Thus, the effective coupling hi appears to be the most important factor in determining 
the end point of the first-order deconfining phase transition line. 

In addition to the transition temperature, other quantities can be estimated using these 
perturbative techniques. For example, the chiral order parameter (ipip) can be estimated 
as a constant term plus a term proportional to the expectation value of a plaquette plus 
a term proportional to the expectation value of the Polyakov loop. However, the use of a 
perturbative evaluation of the fermion determinant obviously fails to include the effects of 
chiral symmetry breaking, which is the dominant factor in determining (ipip) for light quarks. 
Presumably this accounts for the failure of the effective theory for light quark masses. 
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FIGURES 

FIG. 1. Predicted values (circles) of A/3 per fermion at zero temperature versus k compared 
with finite temperature simulation results (squares). 



FIG. 2. Lattice Feynman rules to 0(A 2 ). The Q vertex is not present in the continuum limit. 
FIG. 3. One-loop diagrams contributing to the fermionic renormalization of A/3. 
FIG. 4. Spatial and temporal values of A/3 per fermion versus n for N t = 4, 6, andoo. 
FIG. 5. Typical diagrams from the expansion of A/3 in images. 

FIG. 6. Idealized behaviour of A/3f er mion as a function of /3 pure . 

FIG. 7. Frequency distribution for the expectation value of the Polyakov loop on a 10 3 x 4 
lattice at /3 c (/ti) for hi = 0.025, 0.05, 0.075, and 0.1. 

FIG. 8. A/3 per fermion versus k computed at zero temperature, at finite temperature, and at 
finite temperature with a correction for Polyakov loop effects compared with finite temperature 
simulation results. 
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This figure "figl-5.png" is available in "png" format from: 
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